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Abstract 

This  paper  analyzes  the  problem  of  setting  base-stock  levels  in  a  production  system  that 
operates  under  both  make-to-order  and  make-to-stock  regimes.  We  study  the  problem  of 
choosing  the  base  stock  level  to  satisfy  a  service  level  constraint  and  we  compare  three 
static  priority  schemes  —  first-in-first-out  service,  priority  service  for  make-to-order  jobs, 
and  priority  service  for  make-to-stock  jobs.  Modeling  this  system  by  a  mixed  queueing 
network  and  applying  approximations  derived  from  heavy  traffic  limit  theorems,  we  present 
an  algorithm  for  setting  base-stock  levels  in  each  of  the  three  service  disciplines  and  explore 
conditions  under  which  each  service  discipline  is  favorable. 
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1     Introduction 

The  setting  of  this  paper  is  a  production-inventory  system  that  produces  multiple  types  of 
items.  Some  product  types  are  produced  "to  order";  that  is,  production  of  these  items  are 
initiated  by  actual  customer  orders.  Other  products  are  made  "to  stock";  in  this  case,  the 
system  holds  finished-goods-inventories  from  which  customer  orders  are  filled,  and  production 
of  these  items  is  dictated  by  a  policy  for  keeping  the  inventories  adequately  replenished.  The 
mixture  of  the  two  types  of  production,  make-to-order  (MTO)  and  make-to-stock  (MTS), 
within  a  single  facility  is  representative  of  many  manufacturing  environments,  but  the  inevitable 
compromises  between  capacity  investments,  inventory  costs,  and  customer  service  of  such  an 
operation  are  not  yet  well  understood.  This  paper  explores  two  fundamental  issues:  how  to 
choose  the  base-stock  levels  for  make-to-stock  items  so  as  to  achieve  desired  service  levels,  and 
whether  static  priority  service  policies  can  significantly  improve  the  performance  of  the  system 
in  terms  of  reducing  holding  costs  while  maintaining  satisfactory  service. 

The  stochastic  nature  of  the  production  environment  compounded  with  limited,  finite  pro- 
duction capacity  severely  challenges  the  performance  of  a  production-inventory  system.  Queue- 
ing  network  models,  which  naturally  incorporate  these  two  characteristics,  have  been  employed 
successfully  in  the  study  of  production-inventory  systems  [3,  17].  Past  studies  have  investi- 
gated systems  with  a  large  variety  of  characteristics,  but  to  the  author's  knowledge,  little  has 
been  established  regarding  the  analysis  of  mixed  production  systems  [4,  13,  14].  Building  on 
the  analysis  of  a  queueing  network  model  of  the  mixed  production  system,  this  paper  presents 
guidelines  for  setting  base-stock  levels  as  well  as  priority  schemes  in  such  a  setting. 

We  assume  that  the  production-inventory  system  operates  as  follows.  First,  production  of 
each  type  of  make-to-stock  items  follows  the  "one-for-one  replenishment"  policy.  Under  this 
policy,  a  base-stock  level  is  specified  for  each  product  type;  demands  are  filled  from  finished- 
goods-inventory;  and  each  item  pulled  from  inventory  triggers  a  replenishment  order  to  restore 
the  finished-goods-inventory  to  the  desired  base-stock  level.  Second,  demands  that  cannot  be 
met  due  to  insufficient  inventory  are  considered  to  be  lost.  We  will  investigate  the  performance 
of  the  system  under  three  types  of  static  priority  service  disciplines:  first-in-first-out  service, 
priority  service  for  MTO  products,  and  priority  service  for  MTS  products. 

The  paper  is  organized  as  follows.  Section  2  describes  the  production-inventory  system 
under  study  and  the  corresponding  mixed  queueing  network  model.  In  Section  3  we  present 
formulas  for  setting  base-stock  levels.  In  this  section  we  also  investigate  the  trade-offs  between 
performance  measures  such  as  inventory  level  and  lead  time  for  each  of  the  three  priority 
schemes.  To  simplify  the  presentation,  we  consider  only  production  systems  with  one  make- 
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Figure  1:  A  mixed  production-inventory  system  with  multiple  job  types 

to-stock  job  type  in  this  section.  The  problem  of  setting  base-stock  levels  in  a  production 
system  with  multiple  make-to-stock  job  types  involves  additional  technical  complications  and 
we  return  to  this  problem  in  Section  5  after  we  present  the  theoretical  underpinnings  of  our 
approximations  in  Section  4. 


2     The  Mixed  Network  Model 

The  production-inventory  system  under  study  is  depicted  in  Figure  1.  We  envision  the  produc- 
tion process  as  a  single  operation;  one  can  think  of  the  single  operation  as  an  aggregation  of 
the  entire  production  process;  alternatively,  such  a  characterization  would  be  appropriate  if  the 
production  process  is  consistently  limited  by  a  single  bottleneck  workstation.  The  workstation 
produces  multiple  types  of  products;  make- to-stock  products  are  labeled  as  types  1, . . . ,  d  and 
make-to-order  products  are  given  labels  d  +  1 , . . . ,  c,  so  that  c  is  the  total  number  of  product 
types  in  the  system.  As  stated  in  Section  1,  a  separate  finished-goods-inventory  is  held  for 
each  product  type;  each  make-to-stock  item  follows  the  "one-for-one  replenishment"  policy; 
and  demands  that  occur  when  the  inventory  is  empty  are  considered  to  be  lost. 

We  model  the  make-to-order/make- to-stock  production  system  of  Figure  1  by  the  mixed 
queueing  network  model  with  d  +  1  stations  depicted  in  Figure  2.  Each  station  consists  of  a 
single  server.  Station  0  represents  the  workstation:  an  arrival  at  station  0  signals  a  production 
request  and  each  service  completion  at  station  0  corresponds  to  the  production  of  an  item. 


Class  d+1 


Class  d+2 


Class  c 


Make-to-order 
products 


Make-to-stock 
products 


Class  1 


Class  2 


-Q 


€i 


C/2 


Class  d 


O 


Figure  2:  A  multiclass  mixed  queueing  network 


Stations  1  to  d  model  the  finished-goods-inventories  (FGI)  for  make-to-stock  products:  items 
in  queue  k  represent  FGI  of  type  A;  (1  <  A:  <  d)  and  service  durations  at  station  k  correspond  to 
intervals  between  demands  (i.e.,  inter-demand  times)  of  product  k.  Each  filled  demand  triggers 
a  corresponding  replenishment  order,  so  jobs  that  "depart"  from  station  A;  are  routed  to  station 
0.  Because  demands  that  cannot  be  filled  from  inventory  are  simply  lost,  the  number  of  items 
in  FGI  remains  nonnegative;  moreover,  the  number  of  items  in  FGI  summed  with  the  number 
of  replenishment  orders  at  the  workstation  for  each  product  type  is  constant  at  all  times  and 
equals  the  pre-specified  base  stock  level.  In  the  language  of  queueing  networks,  make-to-order 
products  are  "open"  jobs  whereas  make-to-stock  products  are  "closed"  jobs. 

Demands  for  type  k  products  occur  at  rate  Afc  with  squared  coeflRcient  of  variation  (SCV) 
c1;  their  processing  times  have  mean  rnok  and  SCV  (^q^..  We  do  not  impose  any  restrictions  on 
the  distribution  of  inter-demand  times  and  processing  times;  moreover,  each  product  type  may 
have  its  own  demand  pattern  and  processing  time  distribution.  For  concreteness,  we  will  say 
that  inter-demand  times  and  processing  times  for  each  job  type  form  independent  sequences 
of  independent  and  identically  distributed  random  variables.  This  assumption  is  made  only 
for  ease  of  exposition;  as  our  results  are  based  on  heavy-traffic  limits  of  the  network,  they  are 
valid  under  weaker  conditions  as  well  (see,  for  example,  [7,  9,  10]).  Lastly,  we  will  assume  that 
no  setup  times  are  incurred  when  switching  between  job  types. 

Figure  2  suggests  that  we  model  the  demand  processes  of  make-to-stock  products  by  the 


service  processes  at  stations  1  to  d.  Strictly  speaking,  a  service  process  characterized  by 
independent  and  identically  distributed  service  times  may  not  be  a  faithful  representation 
of  the  demand  process  because  the  first  inter-demand  time  following  a  period  of  inventory 
depletion  typically  is  not  statistically  similar  to  other  inter-demand  intervals.  That  is,  the 
first  service  time  of  a  busy  period  should  be  characterized  by  an  "excess  life"  distribution,  and 
this  is  the  same  as  the  inter-demand  time  distribution  if  and  only  if  the  demand  process  is 
Poisson.  Nonetheless,  the  difference  is  not  significant  in  the  sense  that  the  two  systems  closely 
approximate  each  other  under  heavy  traffic  conditions  (see  [11,  13,  14j). 

We  now  introduce  notation  that  will  be  used  in  the  remainder  of  this  paper.  Write  n^, 
k  =  l,...,d,  to  mean  the  the  target  base-stock  level  of  type  k  make-to-stock  products.  We 
will  assume  that  initially,  the  finished-goods  inventory  of  each  makc-lo-stock  product  type  is 
full.  The  relative  traffic  intensity  at  the  workstation  is  given  by 

c 

Po  =  XI  ^kmok-  (1) 

Let  us  write  pa  (respectively,  pr)  to  mean  the  contribution  from  make-to-stock  (respectively, 
make-to-order)  products  to  the  relative  traffic  intensity  at  the  workstation: 

d  c 

Ps^Yl  ^kTnok  Pr=J2     ^k'^Ok-  (2) 

k=\  k=d+\ 

We  are  interested  in  networks  that  are  neither  pure  open  nor  pure  closed  networks,  which 
implies  Ps  >  0  and  pr  >  0.  Next,  define 

d 

n  =  y]nk     and     Pk  =  —-  (3) 

The  fill  rate  of  type  k  jobs,  defined  to  be  the  fraction  of  type  k  orders  that  are  filled  from 
inventory,  is  denoted  by  Q^t-  The  number  of  type  k  jobs  in  the  workstation  at  time  (  is  denoted 
by  Qok{t),  which  one  can  interpret  as  type  k  work-in-process  (WIP)  inventory.  The  number 
of  items  in  type  k  finished-goods  inventory.  A;  =  l,...,d,  is  denoted  by  Qk{i)  and  satisfies 
the  relationship  Qk{t)  =  n^  -  (5o/t(0-  Finally,  W{t)  is  the  amoimt  of  work  present  at  the 
workstation  at  time  t,  which  includes  the  remaining  processing  time  of  any  job  in  service  at 
that  time. 

Product-Form  Solutions 

It  is  possible  to  derive  explicit  product-form  steady-state  distributions  for  this  network  under 
the  assumptions  that 


(a)  each  station  operates  under  the  FIFO  service  discipline; 

(b)  all  processing  and  inter-demand  times  are  exponentially  distributed;  and 

(c)  all  product  types  have  the  same  mean  processing  time  (at  the  workstation). 
(See  Kelly  [12].)  Under  these  conditions,  it  is  known  that  the  condition 

Pr  <  1  (4) 

is  necessary  and  sufficient  for  stability,  where  pr  is  the  contribution  from  make-to-order  jobs 
to  the  relative  traffic  intensity  at  the  workstation.  Throughout  this  paper  we  will  assume  that 
condition  (4)  holds.  This  of  course  implies  A/t  <  tuq^  for  each  make-to-order  product,  and  we 
will  also  assume  the  same  holds  for  make-to-stock  products.  Moreover,  we  will  assume  that 
make-to-stock  products  do  not  require  full  utilization  of  production  capacity  —  namely,  Ps  <  1 
as  well  —  although  this  is  not  strictly  required  for  stability. 

3     Setting  the  Base-Stock  Levels 

Let  us  begin  by  addressing  the  problem  of  setting  base-stock  levels  in  a  system  with  one  make- 
to-stock  job  type  (d  =  1),  deferring  the  general  discussion  of  multiple  make- to-stock  job  types 
to  Section  5.  (The  number  of  make-lo-order  job  types  is  not  restricted,  so  c  >  2.)  Although 
the  formulas  become  more  complicated  with  multiple  make-to-stock  job  types,  the  behavior  of 
the  system  is  qualitatively  similar. 

Denote  by  ai  the  desired  fill  rate  for  MTS  jobs.  The  problem  is  to  determine  A'^i(Qri), 
N\{ai)  and  Nl{oc\)  the  minimum  base-stock  level  required  to  achieve  this  service  level  under 
FIFO  service,  priority  service  for  make-to-order  jobs  and  priority  service  for  make-to-stock 
jobs,  respectively. 

Base-Stock  Approximations 
Define 

c 

a''  =  Y^  Xkmlkicl  +  c^fc),  (5) 

k=\ 

with  the  interpretation  that  o'^  is  a  measure  of  the  variability  of  the  workcenter's  aggregate 

incoming  workload.  Similarly,  define 

d 
a?  =  5^Afcm^fc(c^  +  c^fc),  (6) 

k=\ 


interpreting  a^  as  a  measure  of  the  variability  of  the  workcenter's  incoming  make-to-stock 
workload.  Let  [xj  denote  the  integer  part  of  a  real  number  x.  We  propose  the  following 
approximations  for  the  base-stock  levels: 


Approximation  1  (Base-stock  Levels) 

I  ^ 

Niiai)    = 


^r(ai) 


J 


Ll-A 


■Niiai) 


N',{ar)     = 


A,a? 


2psil-Ps) 


+  1; 


ai   In     1  + 


1  -  Ps 


Aimoi(l  -  Qi) 


+  1. 


(7) 


(8) 


(9) 


The  formulas  given  in  (7)-(9)  are  based  on  heavy-traffic  analysis  of  the  queueing  network 
model.  This  analysis,  which  we  present  in  next  section,  gives  exact  characterizations  of  the 
system's  dynamics  when  the  workstation  is  balanced  {po  —  1)  and  the  aggregate  base-stock 
level  n  — >  DO,  a  state  of  affairs  that  we  refer  to  as  the  "heavy  traffic  limit."  We  then  modify  the 
exact  formulas  to  obtain  approximations  for  systems  not  in  the  heavy  traffic  limit  with  n  <  oo 
and  Po  not  necessarily  equal  to  one.  The  accuracy  of  these  heavy-traffic  based  approximations 
were  investigated  in  Nguyen  [13,  14],  which  showed  that  for  FIFO  systems,  the  method  provides 
good  approximations  for  performance  measures  such  as  throughput  rates,  inventory  levels  and 
lead  times.  The  next  section  of  this  paper  will  present  further  simulation  experiments,  including 
networks  with  priority  service,  that  also  testify  to  the  accuracy  of  the  approximation  method. 

Encouraged  by  these  preliminary  experiments,  we  now  luni  our  attention  to  the  qualitative 
implications  of  equations  (7)-(9).  We  are  primarily  interested  in  the  relationships  between 
base-stock  levels,  service  levels,  inventory  levels  and  priority  schemes. 

Let  us  first  consider  equation  (7),  the  base-stock  level  under  FIFO  service.  First,  the  base- 
stock  level  is  approximately  proportional  to  the  fill-rate  function  -^^^ ,  which  is  directly  evident 
in  the  case  of  po  =  1  and  can  also  be  verified  for  po  7^  1  by  expanding  the  log  term.  Second,  a 
higher  base-stock  level  is  required  for  a  system  with  higher  aggregate  variability,  as  measured 
by  the  unitless  term 

c 

Third,  the  base-stock  level  increases  when  type  1  utilization,  as  measured  by  Aimoi,  decreases, 
or  equivalently,  because  po  is  approximately  one,  when  there  is  more  competing  work. 
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Taking  FIFO  as  the  "standard"  scenario,  let  us  now  compare  the  base-stock  levels  under 
the  two  priority  disciplines.  Our  approximation  method  reveals  a  very  simple  relationship 
between  FIFO  and  MTO  priority:  To  attain  the  same  service  level  when  priority  is  given  to 
MTO  jobs,  set  the  base-stock  level  to  be  -p^x:  times  that  of  FIFO. 

As  an  example,  suppose  that  MTO  jobs  do  not  command  a  high  portion  of  the  workcenter's 
capacity  (say  1—pr  =  0.90).  If  the  workstation  were  to  give  priority  to  MTO  products,  the  same 
MTS  service  level  can  be  attained  with  little  increase  in  the  base-stock  level  (the  base-stock 
level  would  increase  by  11%),  while  MTO  lead  time  would  be  significantly  reduced. 

When  MTS  jobs  receive  priority,  the  base-stock  level  is  set  as  if  there  were  no  make-to-order 
jobs  in  the  system.  Hence,  the  "aggregate"  variability  a"^  is  replaced  by  its  MTS  component 
<T^,  and  utilization  po  of  the  workstation  is  reduced  to  consist  only  of  MTS  work  pg. 

TYade-offs  Between  Throughput  and  Inventory 

To  get  a  better  a  sense  of  the  trade-offs  between  throughput  and  inventory,  let  us  consider  a 
network  with  one  make-to-order  and  one  make-to-stock  job  type  (d  —  1  and  c  =  2).  We  will 
consider  three  such  systems.  In  all  three  systems,  all  inter-demand  times  and  processing  times 
in  the  network  are  exponentially  distributed;  the  mean  processing  time  for  both  job  types  is 
0.0008;  and  the  utilization  of  the  workstation  is  0.96.  In  the  first  system,  we  set  Ai  =  1000 
and  A2  =  200;  in  the  second  system  Ai  =  A2  =  600;  and  in  the  third  system,  Ai  =  400  and 
A2  =  800.  With  these  parameters,  MTS  utilization  varies  from  0.80  in  System  1  to  0.50  in 
System  2  to  0.32  in  System  3. 

Using  the  approximations  developed  in  Section  4,  we  can  quickly  and  effectively  estimate  the 
performance  of  systems  with  different  parameters  operating  under  different  service  disciplines. 
For  MTS  jobs,  the  primary  performance  measures  of  interest  are  fill  rate  and  average  finished 
goods  inventory.  For  MTO  Jobs,  we  arc  interested  in  the  average  lead  time,  by  which  we  mean 
the  sum  of  the  time  that  an  order  spends  waiting  to  be  processed  and  its  processing  time. 
Rather  than  reporting  the  average  lead  time,  however,  we  will  present  our  results  in  terms  of 
the  "waiting  factor,"  defined  as  the  ratio  of  the  average  waiting  time  to  the  average  service 
time,  which  we  feel  is  a  more  meaningful  characterization  of  the  penalty  due  to  congestion. 

Figure  3  shows  the  average  FGI  and  the  average  waiting  factor  for  each  priority  scheme 
and  each  system  when  a  fill  rate  of  0.90  is  required  for  MTS  jobs.  Each  line  in  the  graph 
corresponds  to  one  of  the  three  systems;  the  left-most  data  point  of  each  line  corresponds  to 
the  poUcy  that  gives  priority  to  MTO  jobs;  the  middle  point  corresponds  to  FIFO  service;  and 
the  right-most  point  arises  from  the  policy  that  gives  priority  to  MTS  jobs.  Figure  4  shows 
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Figure  3:  Attaining  a  fill  rale  of  0.90 

the  analogous  performance  measures  associated  with  a  fill  rate  of  0.95. 

Let  us  investigate  Figure  4  in  more  detail,  starting  with  the  FIFO  data  points.  System  3 
has  the  lowest  percentage  of  MTS  work  so  we  expect  it  to  require  the  lowest  base-stock  for  a 
given  service  level  constraint.  As  seen  in  Figure  4,  the  average  FGI  is  lowest  for  System  3  and 
highest  for  System  1,  which  has  the  largest  percentage  of  MTS  work.  On  the  other  hand,  the 
waiiting  penalty  for  MTO  jobs  is  highest  in  System  3  and  lowest  in  System  1 . 

The  waiting  penalty  for  MTO  jobs  can  be  reduced  by  giving  priority  to  make-to-order  jobs, 
and  the  MTO  data  points  in  Figure  4  show  the  corresponding  increase  in  average  FGI  that  is 
required  to  maintain  a  fill  rate  of  0.95.  The  increase  is  much  more  dramatic  in  System  3  than 
in  System  1,  and  in  fact,  under  MTO  priority,  System  3  requires  the  highest  base-slock  level. 
A  similar  analysis  applies  to  the  MTS  data  points.  If  it  is  expensive  to  hold  FGI  for  MTS 
jobs,  one  can  decrease  the  required  base-stock  level  while  maintaining  the  same  service  level  by 
giving  service  priority  to  MTS  jobs.  This  results  in  a  degradation  of  service  for  MTO  jobs  as 
shown  in  Figure  4,  and  the  degree  of  degradation  in  performance  depends  on  the  parameters 
of  the  particular  system.  Depending  on  the  relative  costs  of  holding  FGI,  costs  of  having  WIP, 
and  costs  of  customer  delay,  it  may  be  desirable  to  operate  System  1  by  giving  priority  service 
to  MTO  jobs,  to  operate  System  2  with  FIFO  service,  and  to  operate  System  3  by  giving 
priority  service  to  MTS  jobs. 

To  get  a  better  appreciation  of  the  impact  of  the  fill  rate  constraint.  Figure  5  shows  the 
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Figure  4:  Attaining  fill  rate  =  0.95 
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Figure  5:  Performance  measure  trade-offs:  System  1 
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Figure  6:  Performance  meaisure  trade-offs:  System  3 

average  FGI  and  average  waiting  factor  in  System  1  for  varying  fill  rates,  ranging  from  0.90  to 
0.99.  As  in  the  previous  graphs,  the  left  most,  middle,  and  right  most  data  points  correspond 
to  MTO  priority,  FIFO,  and  MTS  priority,  respectively.  Figure  6  shows  the  analogous  results 
for  System  3. 

4     Heavy  Traffic  Limits  and  Approximations 

Technical  Preliminaries 

In  order  to  state  the  heavy-trafRc  limit  theorems,  we  need  to  refer  to  a  "sequence  of  systems 
in  heavy  traffic."  Indexing  systems  in  the  sequence  by  n,  we  will  denote  parameters  of  the  n^^ 


,(n) 


(n). 


system  by  a  superscript  (n):  for  example,  A^  '  is  the  demand  rate  of  type  k  jobs  and  Qofc  (0  "s 
the  type  k  work-in-process  inventory  at  time  t  in  the  n'*"  system.  Fixing  a  vector  {0\, . . .  ,0d) 
where  0  <  ^3^  <  1  and  Yfk=\  /5fc  =  1,  we  set  the  base-stock  level  for  type  k  jobs,  k=  1, ...  ,d  in 
the  n'*"  system  to  be  n^  =  /?/cn.  More  precisely,  we  ought  to  choose  an  integer  value,  such  as 
the  integer  part  of  /3fcn  plus  one.  (We  add  one  so  that  we  are  guaranteed  a  positive  base-stock 
level.)  In  the  interest  of  keeping  the  exposition  simple,  however,  let  us  proceed  with  the  former 
prescription  while  keeping  in  mind  that  this  is  without  any  loss  representation  power. 

For  the  limit  theorems,  we  require  that  A^"    — »  Afc  and  tjIq^  — >  mot  as  n  — »  oo.  Moreover, 
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we  will  assume  that  the  following  heavy  traffic  condition  holds: 

n{pQ    -  1)   — »  /i   as  n  — »  oo  (10) 

where  /i  is  a  finite  (but  not  necessarily  nonpositive)  number.  This  condition  implies  that  the 
approximations  derived  from  the  limit  theorems  are  particularly  relevant  when  the  load  at  the 
production  faciUty  is  approximately  balanced  and  the  base-stock  levels  are  large. 

The  heavy  traffic  limit  theorems  will  establish  convergence  of  scaled  processes  as  the  base- 
stock  level  becomes  large,  i.e.,  as  n  — ►  oo.  Specifically,  we  are  interested  in  processes  that  have 
been  scaled  in  the  following  manner: 


iy"(0     =     -W^'^Hn'^t); 
n 

Q^,{t)    =     -Ql^'^^inH),     A;  =  l,...,c; 


n 

Q'iit)  =  l^Qi^\riH),    k  =  i,...,d. 

The  mode  of  convergence  in  the  limit  theorems,  denoted  by  ==>,  is  convergence  of  distributions 
on  the  space  D,  where  D  is  the  space  of  right  continuous  functions  on  [0,  oo)  having  left  limits, 
and  the  topology  is  the  Skorohod  Ji  topology  [2,  16). 

First-in-first-out  service  discipline 

We  will  start  by  reviewing  the  case  of  FIFO  service,  which  was  studied  in  Nguyen  [13,  14]. 
Let  us  assume  that  production  of  all  job  types  follows  the  first-in-first-out  service  discipline. 
Define 

6*  =  min{Ar'/3i,...,A-'/?rf}. 

We  will  write  Qq  and  <5"  to  mean  the  vectors  (Qoi,  •  •  • ,  Qqc)  and  ((57>  •  •  ■ ,  QS),  respectively, 
and  similarly  for  QS  and  Q*. 

Theorem  1  (Nguyen  [14],  FIFO)  Suppose  that  the  production  facility  operates  under  FIFO 
service  and  the  heavy  traffic  condition  (10)  holds,  then  (W^",Qo'Q")  =^  {^*,Qo^Q*)  a^ 
n  — >  oo,  where 

W*  is  reflected  Brownian  motion  on  [0,6*]  with  drift  /i  and  variance  a^,        (11) 
Qokit)  -  ^kW'it),      l<k<c, 
Q'jH)  =  (3k  -  Qo:it),      \<J<d. 
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Remark:     Henceforth  we  will  express  statement  (11)  by  the  abbreviated  notation  "W*  = 

Theorem  1  suggests  that  for  large  base-stock  levels,  the  workload  process  at  the  production 
facility  can  be  approximated  via 

-Win^-)  %  RBM{[0,b'],n{po  -  l),a^). 

Upon  reversing  the  scaling,  we  obtain  the  "heavy  traffic  approximation" 

W{-)^RBMi[0,b],po-l,a''),  (12) 

where 

6  =  minfA^'ni,--- .Aj'rtd}-  (13) 

The  approximation  takes  the  form  of  a  process  that  lives  on  the  bounded  interval  [0,6], 
whereas  the  actual  workload  process  may  assume  arbitrarily  large  values!  This  phenomenon, 
which  was  discussed  in  [13,  14],  may  be  explained  by  saying  that  make-to-stock  jobs,  whose 
volume  is  limited  by  the  base-stock  levels,  effectively  "regulate"  the  total  workload  process  so 
as  to  keep  the  amount  of  work  from  becoming  too  large  relative  to  the  base-stock  levels. 

The  one-dimensional  reflected  Brownian  motion  (RBM)  essentially  captures  all  of  the  dy- 
namics of  the  production-inventory  system,  which  itself  is  c-dimensional  where  c  is  the  number 
of  job  types.  From  Theorem  1,  we  obtain  the  following  approximation  for  type  k  WIP  inven- 
tory: 

Qoki-)^\kW{-),     k  =  l,...,c;  (14) 

and  for  type  k  FGI,  we  have 

Qki-)=nk-Qoki-)    k  =  l,...J. 

Closer  inspection  of  the  above  two  equations  reveals  that  type  A:  PC  1  is  an  RBM  on  the  interval 
[rik  —  Afcb,  rifc).  If  k  is  such  that  6  =  A^'n^,  then  the  approximation  states  that  the  FGI  behaves 
as  a  reflected  Brownian  motion  and  takes  values  between  0  and  n^.  However,  if  /c  is  such  that 
Aj^'njt  >  b,  then  the  approximation  suggests  that  the  FGI  process  is  bounded  away  from  zero. 
In  other  words,  the  heavy  traflric  approximation  suggests  that  type  k  jobs  never  stocks  out! 

Under  the  heavy  traffic  scaling  and  the  eventual  limit,  one  or  more  job  types  will  emerge  as 
"bottleneck  job  types,"  which  are  those  job  types  with  the  smallest  values  of  A^^'ufc.  The  same 
"regulating"  mechanism  that  keeps  the  workload  process  from  becoming  too  large  relative  to 
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the  total  base-stock  levels  also  serves  to  protect  non-bottleneck  job  types  from  ever  stocking 
out.  That  is,  non-bottleneck  FGI's  are  kept  from  becoming  too  small  relative  to  the  FGI's  of 
bottleneck  job  types.  (See  Nguyen  [14]  for  more  details.) 

To  state  our  approximation,  let  us  suppose  that  make- to-stock  product  types  are  numbered 
so  that 

Af'ni  <  Aj 'n2  <  ■••  <  A^'rid  (15) 

(note  that  b  =  A^'ni);  at  the  end  of  this  section  we  will  discuss  the  case  where  some  of  the 
inequalitites  may  be  equalities.  In  addition,  define 


g_2(A)-l) 


(16) 


Approximation  2  (FIFO)  Suppose  that  the  production  facility  processes  jobs  on  a  FIFO 
basis.  The  long-run  average  make-to-order  fill  rates,  make-to-stock  FGI's,  make-to-order  WIP 
inventory  levels,  and  make-to-order  lead  times  are  approximated  by  {a,Q,Qo,T),  respectively, 
defined  as  follows: 

. <72/2moi 


Po  =  1 


ui  -l-cr2/2moi 

1  - f -,,-  ■     I     otherwise. 

{  Aimoi  VI  -e-">«'^/°i'*''/ 


ajt  =  1  - 


1 


P5-1 


fc-1  c 


{k 


,  d),      where 


'^  2(p^  -  1) 

Po  =  ^  OiiXimoi  -\-  ^  A;mo/     and    O''  =       °  ^ — -. 


i=\ 


W  =  < 


;  =  A: 


2Ai  \"'^2moi 
nipo 


Po  =  1 


1 


—  --     otherwise. 


,   aiAi(l  -e-">'^^/°i-^>)      0 

{k  =  l,...,d). 


Po 


Afc, 


Qok  =  —W,      {k  =  d  +  \,...,c). 
po 

fk  =  W-hmok,      (/c  =  d-H,...,c). 


(17) 

(18) 
(19) 

(20) 

(21) 

(22) 
(23) 


Equations  (17)  and  (18)  use  as  data  the  fill  rates  d^,  which  themselves  are  unkown  and 
must  be  approximated.  In  order  for  the  proposed  approximation  to  be  consistent,  we  therefore 
must  show  that  equations  (17)  and  (18)  uniquely  identify  d^,  k  =  l,...,d,  and  that  these 
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parameters  take  values  between  0  and  1.    This  was  done  in  Nguyen  [13].    Moreover,  it  is 
straightforward  to  compute  the  fill  rates  numerically. 

To  end  this  section,  we  return  to  the  question  of  how  to  resolve  possible  equalities  in 
equation  (15).  For  example,  suppose  that  Aj"'ni  =  Aj'n2  <  Aj'ua.  In  such  a  case,  we 
propose  to  approximate  type  2  products  in  the  same  manner  as  type  1  products  (as  if  it  were 
a  bottleneck  type)  and  then  proceed  to  type  3  products  exactly  in  the  same  way  as  outlined. 

Priority  for  make- to- order  products 

We  now  consider  the  case  where  production  gives  priority  to  products  that  are  made-to-order. 
The  following  theorem  is  proved  in  Appendix  A: 

Theorem  2  Suppose  that  make- to- order  products  have  preemptive-resume  priority  in  the  pro- 
duction facility.  Under  the  heavy  traffic  condition  (10),  (M^",<?5,(3")  =>  {W\Ql,Q*)  as 
n  —*  oo,  where 

Q*ok{t)  =  0,       k  =  d+l,...,c, 
W*  =  RBM{{0,Psb*],  /i,  a^), 
Qok(t)  =  —W*{t),     k=l...,d, 
Q;{t)=f3,-Q*o,{t),     j^l...,d. 

Under  the  heavy  traffic  condition  (10)  which  requires  the  production  faciltiy  to  be  balanced 
and  the  base-stock  level  to  be  large,  Theorem  2  states  that  the  (scaled)  number  of  make-to- 
order  products  in  the  network  at  any  time  is  negligible.  Because  make-to-order  products  receive 
preemptive-resume  priority  in  production,  they  are  not  hindered  by  make-to-stock  work  and 
effectively  faces  a  production  center  in  light  traffic.  Therefore  we  find  that  under  the  heavy 
traffic  scaling,  the  number  of  make-to-stock  products  at  station  1  becomes  negligibly  small  in 
the  hmit. 

For  estimation  purposes,  however,  we  propose  that  make-to-order  products  be  approxi- 
mated using  the  standard  heavy  traffic  methodology  for  GI/G/1  queues  that  correct  for  traffic 
intensities  less  than  unity.  Let  us  assume  that  job  are  ordered  so  that  (15)  holds  (with  the 
understanding  that  equalities  would  be  treated  the  same  way  as  described  for  the  FIFO  case). 
Set 

a^^    X:    ^kml,{c^,+cl).  (24) 

k=d+l 

Paralleling  the  arguments  developed  for  the  FIFO  system,  we  arrive  at  the  following  approxi- 
mation: 
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Approximation  3  (Priority  to  Make-to-Order)  Suppose  that  make-to-order  products  re- 
ceive priority  in  production.  The  long-run  average  make-to-order  fill  rates,  make-to-stock 
FGI's,  make-to-order  WIP  inventory  levels,  and  make-to-order  lead  times  are  approximated 
by  {dr,Q'' ,Qq,T^),  respectively,  defined  as  follows : 

a^/2moi 


dT  =  < 


1  - 


(1  -  Pr)ni  +(72/2moi 
I  Aimoi  Vl -e~('"''^)""^^/°i'^» 


Po  =  1 
otherwise 


al  =  l 


Ph" 


AfcTTioifc  Vl  -  e-(*-''")"*''5**'"''/**'^* 


fc-i 


Pq''  =  ^2  Oi'iXimoi  +  Y2  -^;"^oi     and     0 

1=1  l=k 


Tk 


{k  —  2, . . .  ,d),     where 
2(pS'-1) 


W^  =  I 


2Xi{l-pr) 


1 


(1  -  pr)n\ 
(1  -  Pr)poni 


2moi 


{k  =  l,...,d). 


po  =  1 
otherwise. 


Ql  =  n,-'i^W 


{k  =  d  +  l,...,c). 
{k  ^  d+  l,...,c). 


(25) 

(26) 

(27) 

(28) 

(29) 
(30) 
(31) 


In  this  approximation  procedure,  the  make-to-order  products  are  analyzed  as  if  there  are 
no  other  product  types  in  the  system.  That  is,  the  approximations  in  (30)  and  (31)  derive 
simply  from  heavy  traffic  analysis  of  a  multiclass  queue  (see  [9,  10,  15]). 

The  approximations  for  make-to-stock  products  are  based  on  the  results  of  Theorem  2,  and 
the  methods  for  modifying  the  heavy  traffic  limit  results  are  essentially  similar  to  those  of  the 
FIFO  analysis.  Observe  that  in  giving  priority  to  the  make-to-order  products,  the  workload 
process  at  the  production  center  becomes  smaller  (from  Theorem  2,  the  interval  is  now  (0,ps6] 
rather  than  [0,  b]).  This  is  as  we  expect:  by  giving  priority  to  make-to-order  jobs,  make-to-stock 
jobs  are  more  likely  to  be  stocked-out,  the  total  rate  at  which  replenishment  orders  arrive  at 
the  workstation  decreases,  hence  the  workstation  experiences  less  congestion. 

The  relationship  between  queue  length  and  workload  is  also  different  from  the  FIFO  case. 
An  apphcation  of  Little's  Law  to  the  equation  Qok{t)  ^  ^W^(0  suggests  that  we  interpet 
W{t)/ Pa  as  the  waiting  time  of  a  class  k  job  at  the  workstation.  Let  us  replace  p^  by  1  —  pr, 
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noting  that  the  two  are  equivalent  in  the  case  po  =  I.  We  now  have  the  natural  interpretation 
of  1  —  Pr  as  the  amount  of  time  per  unit  time  available  to  the  workstation  for  low- priority  work. 
That  is,  if  a  low-priority  job  arrives  to  see  W{t)  units  of  work  at  the  workstation,  then  it  can 
expect  to  wait,  on  average  W{t)/{1  —  pr)  units  of  time  before  receiving  service. 

Priority  for  make-to-stock  products 

Lastly,  we  consider  the  system  in  which  make-to-stock  products  are  given  priority  in  production. 
The  following  theorem  is  proved  in  Appendix  B: 

Theorem  3  Suppose  that  make-to-stock  products  have  preemptive-resume  priority  m  the  pro- 
duction facility.  Under  the  heavy  traffic  condition  (10),  (M/",<?o.Q")  =>  (t^',<?o.<?*)  «« 
n  — >  CO,  where 

Q'okit)=0,       k  =  l,...,d, 

Ql{t)  =  Pk,       k  =  l,...,d, 

W  ^  RBM{{0,cx)),  po-1,  cr2), 

Qokit)  =  —W*{t),     k  =  d+l...,c. 

POo 

The  limit  theorem  states  that  there  is  no  accumulation  of  WIP  inventory  for  make-to-stock 
products,  which  receive  preemptive  resume  production  priority  and  therefore  experiences  pro- 
duction in  light  traffic.  Moreover,  because  there  is  essentially  no  congestion  in  production, 
finished-goods  inventory  is  always  full  in  the  heavy  traffic  limit.  As  in  the  previous  develop- 
ment, however,  we  propose  that  for  estimation  purposes,  make-to-stock  products  be  modeled  by 
the  standard  heavy  traffic  approximation  for  a  closed  generalized  Jackson  network  [5,  6,  13,  14]. 
In  the  following  approximation,  we  assume  job  types  are  enumerated  so  that  condition  (15) 
holds.  Set 

0^  ,    'AElzA.  (32) 

Approximation  4  (Priority  to  Make-to-Stock)  Suppose  that  make-to-stock  products  re- 
ceive preemptive-resume  priority  in  production.  The  long-run  average  make-to- order  fill  rates, 
make-to-stock  FGI's,  make-to-order  WIP  inventory  levels,  and  make-to-order  lead  times  are 
approximated  by  {a',Q',QQ,f'),  respectively,  defined  as  follows: 

a\  =  i-A-i.    ^:.:^....V  (33) 


X^rnok'  ■        -■"*"*'--"      --      -^--^d.),     where  (34) 
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System 

Ai 

A2 

A3 

moi  -  mo2  =  mo3 

Po 

I 

8.0 

6.25 

2.0 

0.06 

0.975 

II 

10.0 

10.0 

10.0 

0.033 

0.990 

Table  1:  System  Parameters:  Systems  I  and  11 


p''s  =  Yl  "?A/moj  +  Y.  ^'^01     and    9^  =    ^^\ 


d 

z 

l=k 


W^  = 


n\Ps 


afAi(l-e-""''*'/'^J^')      0s' 
Qi^rik-^^W',      {k^l,...,d). 


A) 


Qf 


A/c 


Ok 


2po(l  -  Ps)  V  1  -Po 


{k  -  d+  l,...,c), 


n 


2(l-p,)(l-po) 


+  TUQk,      (A;  =  d  +  1, . . .  ,c),     where 


po  =  Y^afXimoi  +    X!   Ajmoi. 
l=\  l=d+\ 


(35) 
(36) 
(37) 
(38) 
(39) 
(40) 


Note  that  this  approximation  of  make-to-stock  products  is  identical  to  the  approximation 
in  the  FIFO  case  if  we  "eliminate"  all  make-to-order  products.  The  approximation  of  make- 
to-order  products,  which  is  based  on  the  results  of  Theorem  3,  uses  the  same  principles  as 
Approximations  2  and  3.  The  workload  process  now  lives  on  the  entire  non-negative  half-hne. 
Moreover,  rather  than  the  relative  traffic  intensity  po  which  may  be  greater  or  equal  to  one, 
we  propose  to  use  the  actual  traffic  intensity  po  which  will  be  strictly  less  than  one. 

Numerical  Examples 

For  the  case  of  FIFO  networks,  the  approximations  proposed  in  this  section  slightly  differ  from 
those  given  in  our  previous  works  [13,  14).  For  example,  whereas  [13]  approximates  Oi  by 

Aimoi  Vl  -e-"'"/"'^' 
our  present  approximation  contains  an  extra  po  in  the  exponent  term: 

1 L-( ^^^i V 

Aimoi  Vl  -  e-"!^*/**"*"  / 
The  modifications  to  the  approximation  formulas  were  motivated  by  our  heavy  traffic  analysis 
of  priority  service.   Recall  that  the  formulas  presented  in  Section  4  are  asymptotically  exact 
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in  the  heavy  traffic  limit,  which  requires  n  — »  oo  and  Pq  —>  I  Away  from  the  limit,  one  must 
restore  the  finite  base-stock  levels  n*  and  the  utilization  factor  po  to  their  appropriate  roles  in 
order  to  obtain  good  approximations.  For  technical  reasons  that  we  will  not  detail  here,  we 
reahzed  that  the  formulas  are  more  "correct"  when  the  factor  po  is  introduced  in  the  exponent 
term.  This  modification  also  results  in  much  improved  estimates. 

In  addition,  we  proposed  in  this  paper  an  explicit  formula  for  the  fill  rate  of  non-bottleneck 
job  types,  whereas  [14]  established  only  boimds.  To  test  these  approximations,  let  us  consider 
a  network  with  two  MTS  job  types  and  one  MTO  job  type  (d  =  2  and  c  =  3).  We  will  study 
two  systems  whose  service  time  and  inter-demand  time  distributions  are  all  assumed  to  be 
exponential  and  whose  parameters  are  shown  in  Table  1 .  (These  are  the  same  systems  studied 
in  [14].)  Note  that  both  Systems  I  and  11  satisfy  the  product-form  conditions  (see  Section  2) 
so  we  can  compare  our  estimates  against  exact  performance  measures. 

The  performance  measures  of  interest  are  the  throughput  rate  and  FGI  of  make- to-stock 
jobs  and  waiting  time  for  make-to-order  jobs.  Tables  2  and  3  compare  our  approximations 
of  these  performance  measures  against  exact  figures.  In  most  cases,  the  approximations  are 
reasonably  good.  The  method  performs  less  well  when  the  system  is  far  from  the  heavy 
traffic  limit  (for  example,  when  base-stock  levels  are  low)  or  when  there  is  no  clear  bottleneck 
job  type  (i.e.,  when  Af'ni  and  Aj 'n2  are  approximately  equal).  As  one  would  expect,  the 
estimates  improve  markedly  as  the  utilization  of  the  workstation  increases,  as  the  base-stock 
levels  increase,  or  as  a  clear  bottleneck  job  type  emerges. 

To  investigate  the  accuracy  of  the  approximations  for  priority  policies,  let  us  consider  a 
network  consisting  of  two  job  types,  one  make- to-stock  and  the  other  make-to-order  (thus  d  =  1 
and  c  =  2).  We  will  work  with  two  such  systems,  whose  parameters  are  presented  in  Table  4, 
and  for  each  system  we  consider  three  base  stock  levels,  n\  —  10,20, 50.  We  set  the  utilization 
of  the  workstation  to  be  0.90  in  both  systems.  All  processing  times  and  inter-demand  times 
are  exponentially  distributed. 

We  compare  our  estimates  against  those  obtained  from  simulation.  Table  5  displays  the 
results  when  priority  is  given  to  make-to-order  jobs.  Because  of  their  higher  priority,  make-to- 
order  jobs  are  essentially  unaffected  by  the  number  of  make-to-stock  jobs  in  the  network.  (One 
can  think  of  these  jobs  as  experiencing  an  M/M/1  queue.)  The  long-run  average  delay  time 
for  make-to-order  jobs  was  found  to  be  0.13  in  System  III  and  1.63  in  System  IV  for  all  base- 
stock  levels.  Our  approximation  method  did  extremely  well  in  estimating  these  performance 
measures  —  the  estimates  were  virtually  exact  —  as  one  would  expect  from  any  reasonable 
approximation  method,  so  we  did  not  include  these  figures  in  Table  5. 

The  percentages  shown  in  Table  5  have  the  following  interpretations:  the  percentage  as- 
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ni       722 


Fill  Rate:  Type  1  jobs 


HT       Exact     %  En- 


Fill  Rate:Type  2  jobs 


HT 


Exact      %  Err 


System 

1 

5 

10 

85.4% 

85.7% 

-0.4% 

97.9% 

99.3% 

-1.4% 

10 

20 

93.1% 

93.1% 

0.0% 

99.3% 

100.0% 

-0.7% 

15 

30 

95.9% 

95.9% 

0.0% 

99.7% 

100.0% 

-0.3% 

16 

30 

96.2% 

96.2% 

0.0% 

99.6% 

100.0% 

-0.4% 

32 

60 

98.8% 

98.8% 

0.0% 

99.9% 

100.0% 

-0.1% 

48 

90 

99.5% 

99.5% 

0.0% 

100.0% 

100.0% 

0.0% 

16 

20 

96.2% 

96.4% 

-0.2% 

98.8% 

99.7% 

-0.9% 

32 

40 

98.8% 

98.8% 

0.0% 

99.7% 

100.0% 

-0.3% 

48 

60 

99.5% 

99.5% 

0.0% 

99.9% 

100.0% 

-0.1% 

10 

10 

93.1% 

94.4% 

-1.4% 

96.2% 

97.3% 

-1.1% 

20 

20 

97.3% 

97.7% 

-0.4% 

98.5% 

99.2% 

-0.7% 

30 

30 

98.7% 

98.8% 

-0.1% 

99.3% 

99.7% 

-0.5% 

System  ; 

2 

10 

10 

92.2% 

95.2% 

-3.1% 

92.2% 

95.2% 

-3.1% 

20 

20 

96.5% 

97.9% 

-1.4% 

96.5% 

97.9% 

-1.4% 

30 

30 

98.0% 

98.8% 

-0.8% 

98.0% 

98.8% 

-0.8% 

40 

40 

98.7% 

99.3% 

-0.5% 

98.7% 

99.3% 

-0.5% 

50 

50 

99.2% 

99.5% 

-0.4% 

99.2% 

99.5% 

-0.4% 

10 

20 

92.2% 

92.4% 

-0.2% 

98.6% 

99.8% 

-1.1% 

10 

30 

92.2% 

92.2% 

0.0% 

99.6% 

100.0% 

-0.4% 

10 

40 

92.2% 

92.2% 

0.0% 

99.9% 

100.0% 

-0.1% 

10 

50 

92.2% 

92.2% 

0.0% 

100.0% 

100.0% 

0.0% 

20 

50 

96.5% 

96.5% 

0.0% 

99.7% 

100.0% 

-0.3% 

30 

50 

98.0% 

98.0% 

0.0% 

99.6% 

100.0% 

-0.4% 

40 

50 

98.7% 

98.9% 

-0.1% 

99.4% 

99.8% 

-0.4% 

Table  2:  Fill  rate  estimates  for  Systems  I  and  II 
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Til       "2 


MTO  Waiting  Time 


HT      Exact     %  Err 


FGI  at  Station  1 


HT       Exact     %  Err 


FGI  at  Station  2 


HT       Exact     %  Err 


System  1 


5       10 

0.34 

0.33 

2.6% 

2.63 

2.67 

-1.6% 

7.87 

7.60 

3.6% 

10     20 

0.59 

0.59 

0.7% 

5.46 

5.51 

-0.8% 

16.22 

15.93 

1.8% 

15     30 

.83 

.83 

0.0% 

8.51 

8.57 

-0.7% 

24.73 

24.48 

1.0% 

16     30 

.87 

.87 

0.4% 

9.14 

9.20 

-0.6% 

24.45 

24.21 

1.0% 

32     60 

1.44 

1.44 

0.6% 

20.30 

20.40 

-0.5% 

50.75 

50.64 

0.2% 

48     90 

1.82 

1.81 

0.6% 

33.13 

33.24 

-0.3% 

78.33 

78.30 

0.0% 

16     20 

0.87 

0.85 

2.2% 

9.14 

9.26 

-1.3% 

14.50 

14.34 

1.1% 

32     40 

1.44 

1.44 

0.6% 

20.30 

20.41 

-0.5% 

30.78 

30.67 

0.4% 

48     60 

1.82 

1.81 

0.6% 

33.13 

33.24 

-0.3% 

48.34 

48.30 

0.1% 

10     10 

0.59 

0.54 

10.0% 

5.46 

5.73 

-4.6% 

6.34 

6.47 

-2.0% 

20     20 

1.03 

0.99 

4.4% 

11.75 

12.02 

-2.2% 

13.47 

13.54 

-0.5% 

30     30 

1.38 

1.36 

2.1% 

18.80 

19.04 

-1.2% 

21.20 

21.23 

-0.2% 

System  2 


10      10 

0.51 

0.43 

19.0% 

5.27 

5.84 

-9.6% 

5.27 

5.84 

-9.6% 

20     20 

0.92 

0.82 

12.1% 

11.04 

11.84 

-6.8% 

11.04 

11.84 

-6.8% 

30     30 

1.28 

1.18 

9.1% 

17.29 

18.21 

-5.1% 

17.29 

18.21 

-5.1% 

40     40 

1.60 

1.50 

7.2% 

23.99 

24.97 

-3.9% 

23.99 

24.97 

-3.9% 

50     50 

1.88 

1.78 

5.9% 

31.13 

32.11 

-3.0% 

31.13 

32.11 

-3.0% 

10     20 

0.51 

0.51 

-1.4% 

5.27 

5.32 

-0.9% 

14.94 

14.56 

2.6% 

10     30 

0.51 

0.52 

-2.5% 

5.27 

5.30 

-0.5% 

24.90 

24.47 

1.8% 

10     40 

0.51 

0.52 

-2.5% 

5.27 

5.30 

-0.5% 

34.88 

34.47 

1.2% 

10     50 

0.51 

0.52 

-2.5% 

5.27 

5.30 

-0.5% 

44.88 

44.47 

0.9% 

20     50 

0.92 

0.93 

-1.1% 

11.04 

11.09 

-0.4% 

40.74 

40.38 

0.9% 

30     50 

1.28 

1.29 

-0.4% 

17.29 

17.36 

-0.4% 

37.09 

36.78 

0.8% 

40     50 

1.6 

1.59 

1.2% 

23.99 

24.24 

-1.0% 

33.88 

33.88 

0.0% 

Table  3:  Waiting  time  and  inventory  estimates  for  Systems  I  and  II 
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Make-to-stock 

Make-to-order 

System 

A, 

moi 

A2 

mo2 

A) 

III 

1.0 

0.6 

1.0 

0.3 

0.900 

IV 

1.0 

0.2 

1.0 

0.7 

0.900 

Table  4:  System  Parameters:  Systems  III  and  IV 


System  III 

System  IV 

m 

Fill  Rate 

FGI 

Fill  Rate 

FGI 

Sim 

Approx 

Sim 

Approx 

Sim 

Approx 

Sim 

Approx 

10 

0.949 
0.40% 

0.950 
0.11% 

6.31 
1.31% 

6.19 
1.90% 

0.758 
0.83% 

0.532 
29.82% 

5.46 
0.51% 

5.79 
6.04% 

20 

0.990 
0.62% 

0.990 
0.00% 

14.39 

1.02% 

14.19 
1.39% 

0.870 
0.83% 

0.804 
7.59% 

12.57 
0.75% 

12.06 
4.06% 

50 

1.000 

0.67% 

1.000 
0.00% 

42.68 

1.08% 

42.90 
0.52% 

0.965 

0.77% 

0.962 
0.31% 

37.00 
1.16% 

34.93 

5.59% 

Table  5:  Priority  service  for  make-to-order  jobs 

sociated  with  a  simulation  number  corresponds  to  the  95%  confidence  interval  expressed  as  a 
percentage  of  the  nimiber,  whereas  the  percentage  associated  with  an  approximation  is  the  rel- 
ative error  of  the  estimate  as  compared  with  simulation.  The  estimates  of  long-rim  average  fill 
rate  and  FGI  all  fall  within  the  95%  confidence  iterval  for  System  III.  The  approximation  does 
not  perform  as  well  in  System  IV  but  the  estimates  improve  as  the  base  stock  level  increases 
(i.e.,  as  the  system  approaches  heavy  traffic). 

Table  6  compares  approximated  and  simulated  queue  times  for  make-to-order  jobs  when 
make-to-stock  jobs  receive  processing  priority.  As  in  the  previous  scenario,  our  approximation 
performs  extremely  well  in  predicting  throughput  rate  and  FGI  for  make-to-stock  jobs  so  we  do 
not  present  those  numbers  here.  (The  throughput  rates  are  essentially  1.00  for  all  base-stock 
levels;  the  base  stock  levels  are  essentially  n\  —  1.5  and  ni  —  .25  for  System  III  and  System  IV, 
respectively;  and  the  estimates  are  exact  up  to  3  decimal  places.) 

5     Setting  the  Base-Stock  Levels,  Revisited 

The  problem  of  setting  base-stock  levels  is  considerably  more  complex  when  there  are  multiple 
make-to-stock  job  types  because  our  characterization  of  fill  rates  is  not  as  explicit.  In  particular. 
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System  III 

System  IV 

m 

Sim 

Approx 

Sim 

Approx 

10 

10.40 

6.75% 

11.19 
7.60% 

6.46 
4.91% 

6.63 
2.63% 

20 

11.20 
5.35% 

11.25 
0.45% 

6.66 
7.39% 

6.63 

0.45% 

50 

11.40 
6.73% 

11.25 
1.32% 

6.45 
6.50% 

6.63 
2.79% 

Table  6:  Waiting  time  under  priority  service  for  make-to-stock  jobs 


our  estimates  of  fill  rates  depend  not  only  on  the  values  of  the  base  stock  levels  but  also  on 
their  relative  magnitude,  as  measured  by  the  ratios  nk/\k^  where  Uk  and  Xk  are  the  base-stock 
level  and  demand  rate  of  type  k  jobs,  respectively.  These  ratios  determine  "bottleneck"  versus 
"non-bottleneck"  job  types,  and  our  estimation  method  differs  for  a  given  job  type  depending 
on  its  bottleneck/non-bottleneck  status.  This  leaves  us  with  the  unresolved  issue  of  which  type 
to  designate  as  the  bottleneck  job  type. 

Denote  by  a^  be  the  specified  fill  rate  for  job  type  k,  k  =  1,. . .  ,d.  Define  the  ratios 


rk  =  < 


Ctk 


2po(l-Po) 


2Afcmofc  V 1  -  a*: 
ak  In  (  1  4- 


1 


Po 


AfcmofcCl  -  Qk) 


if  Po  =  1 
otherwise. 


(41) 


Let  us  for  now  suppose  that  d  =  2  and  suppose  that  job  types  are  numbered  with  ri  <  r2.  We 
propose  the  following  rule-of-thumb  in  determining  the  bottleneck  job  type: 

(a)  If  Ti  <  0.8r2  then  type  1  is  the  bottleneck  job  type  and  type  2  is  non-bottleneck. 

(b)  Otherwise,  we  consider  ri  to  be  "approximately  equal"  to  r2,  in  which  case  we  think  of 
both  types  as  being  non-bottleneck  job  types. 

The  base-stock  levels  are  then  determined  as  follows: 

(a)  If  A:  is  a  bottleneck  job  type,  then 

Xka' 


Nkiok)    =     ' 


A,a2 


ak 


2po(l-Aj), 


2Xkmok  \\  -ak 
ak  In  (  1  + 


+  1 
-Po 


Afcmofc(l  -Q!jfc) 


if  Po  =  1 


-t-  1    otherwise; 


(42) 
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N:iak)    = 


-I-Pr 


Nkiok) 


+  1; 


(43) 


m<^k)  = 


Xkcr' 


(b)  If  A:  is  a  non-bottleneck  job  type,  then 


ln(^] 


Ok  In    1  + 


l-Ps 


Afcmojfe(l  -Ok) 


+  1. 


(44) 


Nkiak)     = 
N^iak)    = 


1 


Li-A 
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AfcmoA:(l  -ak) 


+  1, 


(45) 
(46) 
(47) 


k  d 

p^  =  ^Q/A;mo/+    ^   A;mo/.        (48) 


l=\ 


Given  the  discussions  in  the  previous  sections,  it  is  straightforward  to  extend  these  approx- 
imations to  many  (d  >  2)  make-to-stock  job  type.  Let  us  therefore  devote  our  attention  to  the 
implications  of  the  base-stock  estimates. 

Table  7  displays  the  base-stock  levels  required  to  achieve  various  service  levels  for  the 
two  systems  described  in  Table  1.  The  desired  service  levels  are  given  in  the  first  column. 
The  second  column  contains  the  computed  base-stock  levels  for  each  job  type,  where  A^;^ 
(respectively,  A^;^)  is  the  base-stock  level  for  type  k  jobs  if  type  k  jobs  were  viewed  as  the 
bottleneck  (respectively,  non-bottleneck)  job  type.  The  third  column  shows  the  computed 
ratios  r^,  and  the  fourth  column  contains  the  recommended  base-stock  levels  using  the  decision 
criteria  outlined  at  the  beginning  at  the  section.  The  exact  base-stock  levels  required  to  achieve 
the  specified  service  levels  were  determined  by  exhaustive  search  (recall  that  these  systems  have 
product-form  solutions)  and  are  shown  in  the  last  column. 

The  designation  of  bottleneck  versus  non-bottleneck  job  type  can  have  significant  impact 
on  the  recommended  base-stock  level.  If  all  jobs  were  treated  as  bottleneck  job  types,  the 
resulting  approximations  would  require  much  higher  base-stock  levels  than  necessary.  On  the 
other  hand,  if  all  jobs  were  treated  as  non-bottleneck  job  types,  the  resulting  recommended 
base-stock  levels  would  not  be  enough  to  deliver  the  required  service  level. 

Moreover,  the  numbers  in  Table  7  illustrate  the  conventional  wisdom  that  it  is  typically 
very  costly  to  require  high  service  levels  for  all  make- to-stock  job  types.    To  be  specific,  let 
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ai,Q2 


Nf     ;vp     N^     N^ 


ri,r2 


Ni,N2     NIN^ 


System  I 


90%, 90% 

8 

6 

7 

5 

0.91,1.10 

6,5 

6,5 

90%,95% 

8 

6 

13 

8 

0.91,1.93 

8,8 

7,7 

90%,99% 

8 

6 

31 

16 

0.91,4.84 

8,16 

7,11 

95%, 90% 

14 

9 

7 

5 

1.63,1.10 

9,7 

9,6 

95%,95% 

14 

10 

13 

9 

1.63,1.93 

10,9 

10,9 

95%, 99% 

14 

11 

31 

20 

1.63,4.84 

14,20 

12,14 

99%, 90% 

35 

20 

7 

6 

4.34,1.10 

20,7 

14,7 

99%,95% 

35 

25 

13 

10 

4.34,1.93 

25,13 

19,11 

99%,99% 

35 

30 

31 

26 

4.34,4.84 

30,26 

28,23 

System  II 


90%,90% 

8 

6 

8 

6 

0.79,0.79 

6,6 

6,6 

90%,95% 

8 

6 

15 

9 

0.79,1.49 

8,9 

6,9 

90%,99% 

8 

6 

46 

20 

0.79,4.55 

8,20 

8,14 

95%,95% 

15 

10 

15 

10 

1.49,1.49 

10,10 

10,10 

95%,99% 

15 

12 

46 

26 

1.49,4.55 

15,26 

14,20 

99%,99% 

46 

36 

46 

36 

4.55,4.55 

36,36 

35,35 

Table  7:  Base-stock  levels 
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us  focus  on  System  II.  The  exact  base-stocks  required  to  achieve  service  levels  of  (90%,  99%), 
(95%, 99%)  and  (99%,  99%)  are  (8,14),  (14,20),  (35,35),  respectively.  (The  recommended 
base-stock  levels  are  similar.)  Note  that  a  higher  service  level  for  type  1  jobs  translates  into 
not  only  a  higher  type  1  base-stock  but  also  higher  base-stock  for  type  2.  In  this  example,  the 
base-stock  level  for  type  2  jobs  experiences  a  more  than  two-fold  increase  when  type  1  service 
level  is  changed  from  90%  to  99%. 

From  the  examples  presented  in  this  paper,  it  is  evident  that  the  problem  of  how  to  set  base- 
stock  levels  in  a  hybrid  make-to-stock/make-to-order  production  environment  involves  complex 
and  subtle  issues.  This  work  provides  simple  and  rehable  estimates  of  performance  measures  for 
such  systems,  allowing  managers  to  quickly  evaluate  the  impact  of  different  operating  scenarios 
and  understand  the  trade-offs  between  service  and  operating  costs. 

A     Proof  of  Theorem  2 

We  will  present  here  an  outline  of  the  proofs  for  Theorems  2  and  3.  The  details  of  the  proofs 
are  typical  of  heavy  traffic  analyses  so  in  the  interest  of  brevity,  we  will  suppress  them  whenever 
possible  and  refer  readers  to  previous  works  for  the  full  proofs.  For  the  purposes  of  the  proofs, 
it  will  be  convenient  to  assiune  that  initially,  the  finished-goods-inventories  are  full  and  the 
workstation  is  empty. 

Because  the  workstation  (station  0)  awards  preemptive  priority  to  open  jobs,  these  jobs 
are  unaffected  by  the  presence  of  closed  jobs  and  experience  the  queue  in  light  traffic.  Lemma 
2  of  Peterson  [15]  applies  directly  to  prove  that  Qofc  ^^  (^  tor  k  =  d+  I,. . .  ,c,  where  C(0  =  0 
for  all  f  >  0.  Equivalently,  denoting  by  Wr{t)  the  amount  of  make-to-order  high  priority  work 
at  station  0  at  time  t,  we  have  W^  =>  C-  We  now  follow  the  same  approach  by  Peterson  [15] 
to  prove  convergence  for  the  processes  of  lower  priority. 

For  fc  =  1, . . . ,  c,  denote  by  Vok{i)  the  sum  of  the  first  i  type  k  service  times  at  the  worksta- 
tion and  by  Aic{t)  the  number  of  class  A;  arrivals  to  the  workstation  by  time  (.  For  k  =  1, . . . ,  d, 
let  Vifc(i)  be  the  simi  of  the  first  i  service  times  at  station  k  following  the  initial  n^  jobs  that 
are  at  station  k  at  time  0.  Set 


E{t)  =  sup 

0<i<l 


E  yok{Ak{s))-s 


interpreting  E{t)  as  the  amount  of  time  available  for  processing  low-priority  make-to-stock 
products  in  the  first  t  units  of  time.  Letting 

E{t)  =  E{t)  -  pst,      and 
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n 
it  follows  from  Peterson  [15]  that  E^  =>  E*  where  E*  is  Brownian  motion  with  drift  /i  (from 
equation  (10))  and  variance  a^  (from  equation  (5)). 

For  k  =  1, . . .  ,d,  let5fc  be  the  renewal  process  associated  service  times  at  station  k;  equiva- 
lently,  Sk  represents  the  inter-demand  process  of  (make-to-stock)  type  k  jobs.  For  k  =  1, . . . ,  c, 
let  5oik  be  the  renewal  process  associated  with  type  k  service  times  at  the  workstation;  for 
ji  =  0, . . . ,  d,  let  Bj{t)  be  the  amount  of  time  station  j  is  busy  in  [0,  t]\  and  let  Ij{t)  =  t  —  Bj{t) 
be  the  cumulative  idleness  process  for  station  j.  If  we  let  Tkit)  be  the  amount  of  time  the 
workstation  has  devoted  to  type  k  service  during  the  first  t  units  of  time,  it  follows  from  the 
previous  definitions  that 

d 

^Tk{t)  =  E{t)-Io(t). 

For  this  section,  let  us  denote  by  Wk  the  workload  process  at  station  k  {k  =  0, . . .  ,d).  The 
amount  of  work  jissociated  with  the  initial  Uk  jobs  at  station  A;  (A:  =  1, . . .  ,d)  is  accordingly 
denoted  as  Wic{0).  We  have 


^d 


(49) 


^,^.  ^  I    ELi  VokiSkiBkim  +  ELd+1  VokiAkit))  -  Bo{t)         if  fc  -  0, 

\   WkiO)  +  Vk{Sok{n{t)))  -  Bk{t)  \^k^l,...,d. 

Let  us  now  define  the  "centered"  processes 

Vok{t)     =    Vok{t)-Tnokt  Sokit)     =    Sok{t)-{l/mok)t, 

Vk{t)    =     Vk{t)-{l/\k)t  Skit)     ^    Sk{t)-Xkt, 

Ak{t)    =    Ak{t)-\kt  tkit)     =    Tk{t)  -  XkTnokt. 

Define  T]{t)  to  be  the  arrival  time  to  the  workstation  of  the  low  priority  job  receiving  service 
at  station  0  at  time  t,  and  set  r]{t)  =  i  if  no  low  priority  job  is  in  service  at  that  time.  The 
amount  of  low-priority  time  available  between  7j(t)  and  t  is  exactly  the  amount  of  low-priority 
work  present  at  time  T]{t)  less  the  remaining  service  time  of  the  low  priority  job  currently  in 
service,  if  there  is  one.  That  is, 

E{t)  -  E{r,{t))  =  [Woivit))  -  Wrivim-eiit), 

where  ei{t)  represents  the  remaining  service  time  if  the  job  currently  in  service  is  low-priority, 
and  is  zero  otherwise.  In  terms  of  centered  processes,  we  have 

Psit  -  v{t))  =  [e(77(0)  -  m]  +  \woivit))  -  Wrivm  -  f  i(t)- 


26 


Similarly,  the  allocation  processes  obey  the  property 

Tk{t)  =  VokiSk{Bkir]it))))  +  €2kit),      l<k<d, 


(50) 


where  62* (0  is  the  amount  of  service  the  current  job  has  received  if  that  job  is  of  type  A:  and 
is  zero  otherwise.  Writing  (50)  in  terms  of  centered  processes,  we  arrive  at 


Tkit) 


Vok{SkiBk{v{tm+mokSk{Bk{r,it)))  -  Xkmoklkirjit))  - 
Afcmofc 


{Eivit))  -  m  +  Wo{ri{t))  -  Wr{r){t))  -  ex{t)]  +  e2k{t).  (51) 


Moreover,  (49)  can  be  written  as 

'  Ef=i  [Vok{Sk{Bk{t)))  +  mokSk{Bk{t))  -  Xkniokhit)] 


Wk{t)  =  <^ 


+ 


VokiAkit))  +  mokAk{t)\  +  n(po  -  l)t  +  Io{t) 
Wk{0)  +  Vk{Sok{Tk{t)))  +  ^Jok{Tk{t))  + 


if  A:  =  0, 


if  /c  =  1, . . .  ,d. 


Let  us  introduce  the  following  scaling  conventions:  for  X^^^  a  generic  process,  then 

X"(t)  =  -X^'^Hnh)       X"(t)  -  -^X(")(n2f). 
n  n^ 

Setting 

ELi  \Vok{Sk{Bk{t)))+mokSk{Bk{i))\  + 


Ut)  =  { 


Efc= 


d+\ 


Vok{Ak{t))  +  mokMi)    +  "(Po  -  l)t 


if /c  -0, 


Vk{Sok{Tk{t)))  +  ^Jok{Tk{t))  + 


•^*"»o* 


^  {Vok{Sk{Bkm))))  +  m^kSk{Bk{Tj{t)))}  - 
^  [E{r^{t))  -  m  4-  Q{r){t))  -  W,{r^{t))  -  tx{t)]  + 


J_ 

■mi 

the  scaled  workload  processes  can  be  written  as 

d 

^(n)^('»)  jn 
k=\ 


if  A:  =  1,. . .  ,d. 


(52) 


(53) 


k=l 

w;:{t)  =  Ckit) - ^  (isit)  +  Y. >^i^^rn^^k^kit)]  +  mt) - /?(^"(i))i ,  k  =  h... m) 

Ps    '     \  k=l  I 

/"(■)  is  continuous  an  nondecreasing  with  /"(O)  =  0,  j  =  0, . . . ,  d;  (55) 

/"(•)  increases  only  at  times  t  when  W"(t)  =  0,  j  =  0, . . . ,  d.  (56) 

The  convergence  of  W^  to  the  desired  limits  follows  from  the  results  of  Nguyen  (14).  The  con- 
vergence of  scaled  queue  length  processes  (^ot  and  Q\  are  proved  similary  with  the  methodology 
of  Peterson  [15].  .     D 
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B     Proof  of  Theorem  3 

Closed  jobs  receiving  preemptive  priority  at  station  0  are  not  affected  by  open  jobs,  so  the 
dynamics  of  closed  jobs  in  this  system  may  be  described  a  a  closed  network  with  type  k  jobs 
circulating  between  station  0  and  station  k.  The  results  from  Nguyen  [14]  can  therefore  be 
applied  directly  to  obtain  convergence  for  the  scaled  queue  length  processes  Qq^  and  Q^, 
k  —  I, . . .  ,d.  Turning  to  the  total  workload  at  station  0  and  using  the  notation  developed  in 
Appendix  A,  we  have 

d 

The  results  from  analysis  of  high  priority  make-to-stock  jobs  imply  I]^  ^=>  ^  where  C(0  =  0 
for  all  A;  =  1, . . . ,  d  and  (  >  0.  We  can  therefore  can  conclude  that  W^  converges  to  an  RBM 
on  the  half  line  (0,  oo)  with  drift  fi  and  variance  a"^.  Convergence  for  the  queue  length  of 
make-to-order  jobs  follow  from  the  methods  of  Peterson  [15].  □ 


28 


References 

[1]  Baskett,  F.,  Chandy,  K.  M.,  Muntz,  R.  R.,  and  Palacios,  F.  G.  Open,  closed  and  mixed 
networks  of  queues  with  different  classes  of  customers.  Journal  of  the  ACM  22,  248-260 
(1975). 

[2]  Billingsley,  P.  Convergence  of  Probability  Measures.  Wiley,  New  York,  1968. 

[3]  Buzacott,  J. A.  and  Shanthikumar,  J.G.  Stochastic  Models  of  Manufacturing  Systems. 
Prentice  Hall,  Englewood  Cliffs,  1993. 

[4]  Carr,  S.  A.,  Giillii,  A.  R.,  Jackson,  P.  R.  and  Muckstadt,  J.  An  Exact  Analysis  of  the  No 
B/C  Stock  Policy,  preprint,  1993. 

[5]  Chen,  H.  and  Mandelbaum,  A.  Stochastic  discrete  flow  networks:  Diffusion  approximation 
and  bottlenecks.  Annals  of  Probability  19,  1463-1519  (1991). 

[6]  Dai,  J.  G.  and  Harrison  J.  M.  The  QNET  method  for  two-moment  analysis  of  closed 
manufacturing  systems.  Submitted  for  publication  (1992). 

[7]  GljTin,  P.  W.  Diffusion  Approximations.  In  Handbook  on  OR  and  MS,  Vol  2.,  D.  P. 
Heyman  and  M.  J.  Sobel  (eds.),  Elsevier  Science,  North  Holland,  145-198  (1990). 

[8]  Harrison,  J.  M.  Broxvnian  motion  and  stochastic  flow  systems.  Wiley,  New  York,  1985. 

[9]  Harrison,  J.  M.  and  Nguyen,  V.  The  QNET  method  for  two-moment  analysis  of  open 
queueing  networks.  Queueing  Systems:  Theory  and  Applications  6,  1-32  (1990). 

[10]  Harrison,  J.  M.  and  Nguyen,  V.  Brownian  models  of  multiclass  queueing  networks:  Cur- 
rent status  and  open  problems.  Queueing  Systems:  Theory  and  Applications  (forthcom- 
ing). 

[11]  Iglehart,  D.  L.  and  Whitt  W.  Multiple  channel  queues  in  heavy  traffic  I.  Advances  in 
Applied  Probability  2,  150-177  (1970). 

[12]  Kelly,  F.  P.  Reversibility  and  Stochastic  Networks.  Wiley,  New  York,  1979. 

[13]  Nguyen,  V.  Fluid  and  Diffusion  Approximations  of  A  Two-Station  Mixed  Queueing  Net- 
work. Mathematics  of  Operations  Research,  to  appear. 

[14]  Nguyen,  V.  A  Multiclass  Hybrid  Production  Center  in  Heavy  Traffic.  Submitted  for  pub- 
lication. 


II 


MIT  LIBRARIES 


3  9080  00932  7815 


[15]  Peterson,  W.  P.  A  heavy  traffic  limit  theorem  for  networks  of  queues  with  multiple  cus- 
tomer types.  Mathematics  of  Operations  Research  16,  90-118  (1991). 

[16]  Whitt,  W.  Some  useful  functions  for  functional  limit  theorems.  Mathematics  of  Operations 
Research,  5  67-85  (1980). 

[17]  Yao,  D.D.  Ed.  Stochastic  Modeling  and  Analysis  of  Manufacturing  Systems.  Springer- 
Verlag,  New  York,  1994. 


;|97     '!8r 


Date  Due 


fc  '  ■•<  199{ 


Lib-26-67 


